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METHOD  FOR  COMPFNSATING  TRANSPORT  LAGS  IN  COMPUTER 
IMAGE  GENERATION  VISUAL  DISPLAYS  FOR  FLIGHT  SIMULATION 


I.  INTRODUCTION 

One  of  the  major  problem  areas  in  (light  simulation  today  is  integration  of  visual,  motion,  and 
aerodynamic  subsystems.  A particular  subproblem  is  the  compensation  of  the  visual  subsystem  to 
correspond  to  the  aerodynamic  (usually  called  “flight”)  system.  This  compensation  is  normally  both 
software  and  hardware  controlled.  In  Computer  Image  Generation  (C1G)  systems,  a fixed  transport  lag  on 
the  order  of  50  to  100  milliseconds  needs  to  be  compensated.  Alter  compensation  is  achieved,  the  resulting 
output  is  smoothed  to  remove  flutter  effects.  This  paper  introduces  a general  method  for  this 
compensation. 


It.  PROBLEM  EORMl  I.ATION 


I Basic  Assumptions 

a.  rite  primary  iteration  interval  of  both  the  (light  and  visual  systems  is  h (a  typical  value  for  h is 
1/30  second).  Hi  is  is  not  a limitation  to  the  method,  but  certain  modifications  must  be  made  if 
interpolated  output  to  the  visual  processor  is  used. 

b.  After  each  lliglil  output,  position,  velocity,  and  acceleration  data  are  current  for  all  orientation 
and  position  variables. 

c.  The  llight  output  is  considered  accurate:  that  is.  we  treat  the  simulation  mathematical  model  as 
though  it  were  representative  of  the  actual  aircraft. 

d.  The  sequence  of  future  accelerations  and  the  position,  velocity,  and  acceleration  at  any  moment 
in  time  serves  to  define  the  entire  position  and  velocity  history  forward. 

2.  Integration 

From  mu  assumptions,  we  have  the  following  integration  model  for  arbitrary  variable  v (here 
denotes  the  usual  differentiation  with  respect  to  time). 
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Where  (KKl  for  illustration  purposes. 
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Assumption  (Id)  implies  that  the  knowledge  of  a,,  a?,  a together  with  v|(,  vn,  and  v determine 

Vn+i’  vn+2’  vn+2 +t’  Vn+ 1 ’ vif+2’  311  ^ Vn+2  +t'  ^ons'dcr,  for  example,  the  trapezoidal  method  of 
integration,  where: 


v . . = v +h[v  , +v  I ' 

n+  1 n - 1 n+ 1 n 1 

2 


v . . = v -hO  |v  ,+v  I 
n+ 1 n 2 n+ 1 n * 


This  method  is  clearly  of  the  type  specified  by  paragraph  Id. 

Figure  1 highlights  another  important  property  of  CIC.  compensation,  namely,  that  the  final  variable 
to  be  predicted  (in  this  case  vj’+,+T)  may  not  lie  evenly  on  an  integral  multiple  of  h.  For  our  purposes,  it  is 
necessary  to  consider  the  integration  operator  as  a continuous  one.  Normally,  we  ms.  retain  the 
acceleration  derivative  sequence  values  up  to  any  desired  length  and  the  prediction  point  in  the  future  may 
require  extrapolation  forward  any  number  of  iteration  intervals.  However,  many  prediction  techniques  also 
induce  flutter,  the  rapid  oscillatory  response  of  a visual  system  to  high  frequency  pilot  control  inputs: 
therefore,  it  is  desirable  to  couple,  in  line,  a low-pass  filter  to  remove  high  frequency  effects.  This  leads  us 
to  the  formal  definition  of  the  compensation  problem. 


3.  CIG  Transport  Compensation  Definition 

For  each  variable,  v.  required  by  the  visual  system  for  display  purposes. 

Gwen  vn,  vn,  vn,  vn_, ....  v^.,  find  coefficients  c, , cr  c, cs+3  such  that 

Pn=ClV„+Vn+  £ C3+iVi 


is  a “smooth”  approximation  to  vj’+k+T 


where  k,  and  o<t^1  are  given. 

Here  the  only  open  ended  definition  is  that  of  “smooth."  Some  researchers  have  defined  a more 
training  oriented  definition  (Ricard,  Norman,  & Collyer,  lc)76).  However,  the  generali/ability  of  our 
technique  is  partially  due  to  the  capability  of  adjusting  this  smoothness  parameter  for  each  degree  of 
freedom  of  aircraft  motion. 


Problem  Solution 

Since  the  integration  technique  itself  is  completely  deterministic,  we  can  find  constants.  1,1,.  |v  |4. 
. . .It + 4 such  that 


Vn  + k+T  ~ ,lVn  + ,2Vn  + l.l\r  + ,4ai+,Sa2+  ‘ • ■lk  + 3ak  + ,k  + 4 a 

where  the  parameters  aj,  a2, . . .,  ak,  a are  determined  as  linear  combinations  of  the  vf . v ( 
Ve  > an<1  coefficients  I, , I2>  etc.  depending,  in  general,  on  h. 

First, 

as  = vn  +s  VVn+s(s+l)y3vn+  . . . s(s+I ) (s+2)(  . . . ) (s+M ) y V 
2!  • 6! 
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where  s=l,  2,  3,  . . . k,  k+T(ak+T  = a)  ;(Hildebrand,  1956) 

however,  even  though  the  above  relation  can  be  expressed  conveniently  in  matrix  operator  form,  a 

modification  so  that  the  vector  (vn,yvn,  . . . ,ycvn  ) is  transferred  into  (vn,  vn_i,  vn_2 vn.c)  's 

required.  Using  the  notation  of  Hildebrand  (1956)  wherein 


(;)"•(;)-  ■(') 


= s(s+l)etc, 

2 T 


we  can  Write  the  vector  (a, , a.  ....  a,  , a)  as  the  product  of  two  matrices  and  the  vector  (v  , v , , v ,, 
..  I l k ' * vn’n-l  n-2 
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General  Term:  / -t 
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where  tj  = i i=l,  2, . . . , k 


‘k  + l = k+T 


Imbedding  v into  v by 
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and  a into  a by 
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Here  we  imbed  the(k+l)  X(S+1)  product  T*B  into  the  (k+4)  X (if+3)  product  M 
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Simplifying  notation,  we  let 
M = + M2 

implying  that 
■+  -*■ 
a = MV 

n n 
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with  these  conventions, 

v!,'+ur=i'Mvr 

finally . letting  ('  - I 'M 
we  have 


V1’  =f'V 

n + k+7  1 Vn 

where  denotes  transpose.  That  is.  the  scalar  estimate  of  v| +k+T  is  just  the  inner  product  of  a fixed 
vector  C and  the  vector  Vn,  of  held  data.  The  selection  of  C (equivalent  to  a bounded  linear  functional  over 
the  Hilbert  space  for  which  is  a vector)  is  precisely  determined  by  the  topology  induced  on  the  space 
containing  Vn  by  the  integration  method  used.  The  remaining  structure  on  the  space  is  defined  by  the 
low-pass  filtering  requirement.  Although  several  lirst  order  filters  arc  available,  we  have  chosen  Butterworth 
filters.  These  filters  possess  the  advantages  of  smooth  magnitude  and  phase  response.  The  interested  reader 
should  consult  Oppenheim  and  Shale rt  (1975).  and  Stems  (1975)  lor  additional  information  regarding 
digital  filtering  methods. 

A p 

I or  notation al  convenience,  we  adopt  the  notation  I’  = ^ n + i,*7  as  1 demanded  position  prediction 
at  time  nh  lor  the  CK is.  and  p*  as  the  optimum  (filtered)  position.  The  two  arc  related  by 

p*+i  = fp*+iL£L(ptl+i  +p„i 

where  (Xf<l  is  tlie  Butterworth  filler  constant 

Substituting  pn  = C Vn,  pn+ , = CV  i+ , gives 

pn+i  =fP*  + (y)f(vi,) 

with  initial  condition  given  by  p = p = C V 

Due  to  the  simplicity  of  the  above  form,  it  may  easily  be  executed  either  by  software  or  hardware. 

Remarks  and  Recommendations 

With  appropriate  modifications,  this  technique  can  be  extended  to  combinations  of  transport  and 
phase  lag.  as  are  common  in  model-probe  visual  systems  Due  to  its  gcncralizability,  off-line  evaluations  of 
candidate  solutions  to  the  overall  integration  problem,  taking  into  account  differing  iteration  rates, 
servo-characteristics,  and  simulation  requirements  is  possible.  Research  conducted  at  the  Naval  Training  and 
Equipment  Center  (NThC)  has  demonstrated  the  positive  value  of  sueha  technique  (Ricard et  al..  1976).  We 
recommend  the  use  ol  this  analytical  tool  in  the  off-line  analysis  of  full  mission  flight  simulators  for 
maximizing  cue  correlation. 
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